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Abstract 

Stationary and nonstationary, in particular, chaotic states in long Josephson 
junctions are investigated. Bifurcation lines on the parametric bias current- 
external magnetic field plane are calculated. The chaos strip along the bi- 
furcation line is observed. It is shown that transitions between stationary 
states are the transitions from metastable to stable states and that the ther- 
modynamical Gibbs potential of these stable states may be larger than for 
some metastable states. The definition of a dynamical critical magnetic field 

characterizing the stability of the stationary states is given. 
74.50+r, 05.45.+b 
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INTRODUCTION 



Dynamical chaos in long Josephson junctions is of great interest because it can be a 
source of dynamical noise in devices based on them, in particular, in SQUIDs, limiting the 
sensitivity of these devices. Furthermore, dynamical chaos in long Josephson junctions (LJJ) 
is a very interesting physical phenomenon taking place in nonlinear systems in the absent 
of an external stochastic force. Dynamical chaos in a LJJ is easily excited and therefore 
it may also be investigated experimentally rather easily.00 

In our previous works@'@ we have shown that among a set of solutions of the Ferrell- 
Prange equation describing stationary states of the LJJ in an external magnetic fieldS are 
both stable and unstable ones. At the same time, these stationary states are asymptotic 
solutions of the nonstationary sine-Gordon equation and we have also shown that a selection 
of the stable solutions can be governed by a rapid damping in time of the initial perturba- 
tion entering into the nonstationary sine-Gordon equation through the boundary conditions. 
Changing the intensity of this perturbation at fixed shape, we can obtain various station- 
ary states for the LJJ without a bias current or three clusters of states (stationary, and 
time dependent regular and chaotic) in the presence of a bias current. It turned out that 
asymptotic states are very sensitive to an external perturbation, its value and shape define 
the state (stationary, regular or chaotic) to which the system will tend at t — > oo (we have 
called this influence on the selection of asymptotic states of the small rapidly damping initial 
perturbation in time an effect of memory). The fact of coexistence of all these three charac- 
teristic asymptotic states selected only by the form of the initial perturbation seems to be 
astonishing. It is evidently enough that the Ferrell-Prange equation will not have solutions 
at a large bias current j3. Therefore the question arises: at which values of (3 do stationary 
states of a LJJ disappear or what will be a boundary in the parametric (3 — H$ plane {Hq is 
an external magnetic field) that separates this plane on the regions where stationary states 
do and do not exist? Since the number of solutions of the Ferrell-Prange equation changes at 
variation of the parameters (H , (3), another question arises: what is the form of bifurcation 
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lines in the plane (3 — Hq that separate the parametric plane on the regions with a different 
number of stationary solutions of the Ferrell-Prange equation? 

The existence of several stable solutions of the Ferrell-Prange equation is equivalent to 
the fact that thermodynamical Gibbs potential G associated with the distribution of the 
magnetic field along the junction has minima, and each minimum corresponds to a certain 
solution of the Ferrell-Prange equation. Does a global minimum of G correspond to the most 
stable state (e.g., in the Lyapunov sense)? In the case of the junction of the finite length both 
Meissner and one-fluxon states are thermodynamically advantageous simultaneously, so it 
is interesting to investigate dynamical properties of these states. Answering this question, 
we introduce a dynamiclal critical field that describes the stability characteristic of the 
junctions. 

In Sec. 1 bifurcation lines on the parametric (3 — H Q plane are calculated. In Sec. 2 the 
definition of the dynamical critical magnetic field is given and the dependence of this field 
on /3 and the length of the junction L is calculated. In Sec. 3 transitions between states are 
described. It is shown in Sec. 4 that a chaos strip arises along the bifurcation line on the 
parametric (3 — H plane. The last Sec. 5 contains the discussion of our calculation and brief 
conclusions. 

I. BIFURCATION LINES 

Stationary states of a LJJ are investigated using the numerical integration of the Ferrell- 
Prange equation: 

(Pxx(x) = sin<p(x) - P, (1) 

where (p(x) is the stationary Josephson phase variable, (3 is the dc bias current density 
normalized to the critical current j c , x is the distance along the junction normalized to the 
Josephson penetration length Aj = ^JC^o/8n 2 j c d, $ is the flux quantum, d = 2Xl + b, Xl 
is the London penetration length, b is the thickness of the dielectric barrier. The boundary 
conditions for Eq. (P have the form 
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Vx{x)\x=0 = <Px{x)\x=L = -Hq, (2) 

where L is the total length of the junction normalized to A j and H is the external magnetic 
field perpendicular to the junction and normalized to H = $o/2vrAj<i. 

Numerical integration of Eqs. (|]])— (0) allows us to find the regions with a certain number 
of solutions on the parametric (3 — H plane (Fig. |ip. It is easy to show that the set of 
points corresponding to the even number of solutions forms two-dimensional domains on 
this plane, whereas the set corresponding to the odd ones may form just one-dimensional 
curves. Mostly, the lines corresponding to the odd number of the solutions of the Ferrell- 
Prange boundary problem coincide with the bifurcation lines. Using the shooting method 
for solving of the boundary problem one can prove that the 27r-periodicity of the function 
H(ip ) expressing the dependence of the magnetic field at the right side of the junction 
(x = L) on the phase taken at the left side (x = 0) results in the appearing of the odd 
number of solutions only when the H(ip ) touches the line H = H in an extreme point, 
i.e., dH(<p )/d<p = 0. As an illustartion, we have plotted in Fig. ^ the function H({p ) at 
H = 0.5, L = 5, p = 0.25 and (3 = 0.45. 

Boundaries between the regions - bifurcation lines - define an essential modification of 
the system. The bifurcation lines in Fig.l are obtained for L = 5; here a step by (3 is equal 
to 5 • 10~ 3 and a step by Hq is equal to 2.5 • 10~ 3 . In this figure the numbers of solutions 
of Eq. (0)-© are pointed out, the numbers of stable solutions are given in the brackets, 
while M and If denote a stable Meissner and one-fluxon states, respectively. It is seen 
that a Meissner state is stable at small values of H and at large values of H a one fluxon 
state is stable. It should be noted that the region where there are no stationary solutions 
(region 0) bounds with the region having a minimum of stationary solutions, being equal 
to 2 (region 2). In approaching the boundary of region and 2 the number of stationary 
solutions decreases: 6 — > 4 — » 2 — >• 0, on the other hand, a number of nonstationary states 
which are the asymptotic solutions of the sine-Gordon equation, increases. Our calculations 
have shown that one of two stationary solutions in region 2 is stable, and another is unstable 
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(metastable). We noted earlieilli that the stable states are symmetrical. The presence of 
bias current (3 leads to a symmetry violation that results, evidently, in the instability of the 
states. 

The problem of the stability of stationary states (p(x) was solved in the following way:0 
the sine-Gordon equation was linearized in the vicinity of stationary solution: <p(x,t) = 
f(x) + 9(x,t), where 9(x,t) is the infinitesimal perturbation. The equation for 9(x,t) - the 
linearized sine-Gordon equation - we can solve by means of the expansion of this function 
in terms of a complete system of eigenf unctions of the Schrodinger operator with potential 
cos [<£>(x)]: 

6(x,t) = J2e Xnt u n (x), (3) 

n 

where u n (x) are eigenfunctions of the Schrodinger operator of the problem: 

— u xx (x) + u(x) cos^(x) = Eu(x), (4) 

Ux{x)\ x =Q = Ux{x)\ x =L = 0, 

and 



A n = - 7 ± vV - E n , (5) 

where 7 is the dissipative coefficient in the sine-Gordon equation. We note that values of A 
coincide with corresponding values of Lyapunov exponents in the case when perturbations 
are considered with respect to the stationary solutions. In general case, Lyapunov exponents 
are calculated in the same way as in Ref. 13. Thus, in the presence of a bias current we 
have the different picture of a LJJ states than at (3 = (this case has been examined in Ref. 
12). For example, at H = 1.9 the increasing of /3 from to 0.22 leads to the changing of 
the stationary states number 6 — > 4 —>■ 2 — >• 0, i.e., to a consecutive losing of the stationary 
solutions. Simultaneously, an increasing of the number of nonstationary states occurs that 
we found by directly solving the nonstationary sine-Gordon equation. 



II. DYNAMIC CRITICAL FIELD 



In the literature the critical magnetic field H Cl in a LJJ is defined as a field value, at 
which an existence of a Josephson vortex (fluxon, soliton) becomes advantageous thermody- 
namically for the first time (see, for example, Refs. 10 and 11). In the case of an infinitely 
long junction the critical field is H cl (oo) = A/tt ~ 1.274. Essentially, this field corresponds to 
the global minimum of the thermodynamic Gibbs potential for the one-fluxon state. How- 
ever, in a junction of finite length there are some local minima that coexist with the global 
one and every minimum corresponds to the solution of Eqs. (|T])-(|2]). Some of these solutions 
are stable, another unstable in the sense discussed in Sec. 1. 

We write down the thermodynamic Gibbs potential in the form 

r L 1 

G= dx[-ip x 2 (x) + 1 - cos^(x) - 0ip(x) - H ip x (x)]. (6) 
Jo I 

Here G is the thermodynamic Gibbs potential per unit length along an external magnetic 
field and normalized to G — $o/167r 3 A jd. The Ferrell-Prange equation is an extremal of 
the functional (|6]). An investigation of the second variation of G shows that all extrema of 
this functional satisfy to the necessary and sufficient conditions of strong miuimum.i Thus, 
all solutions of Eqs. (|I])-(Q) (both stable and unstable ones) correspond to minima of the 
thermodynamic Gibbs potential; one of them is global, the others are local. Our calculations 
of the thermodynamic Gibbs potential (P) show that, for example, at (3 = 0, L = 5 and 
H = 0.67 the Meissner state has a global minimum (Gm = —0.44), but the stable one- 
fluxon state has a local one (Gi/ = 4.03). The one-fluxon state has a global minimum of G 
starting at H = 1.57 (G\f = —2.582) and at the same value of (3 and L. At this value of 
Hq a Miessner state has a local minimum Gm = —2.58. At Hq > 2.09 the Meissner state 
disappears. Thus, at a field less than the critical one H c i, the stable one-fluxon state exists. 
We shall further call the minimum value of a magnetic field at given L and /3, at which 
the stable one-fluxon state appeares for the first time and which corresponds to the local 
minimum of the thermodynamics Gibbs potential as the dynamical critical field H^. It is 
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interesting that the dynamical critical field Hd c makes up on the parametric plane a line 
that coincides with the bifurcation line BC (see Fig. [I]). Our calculations show that the bias 
current increases the dynamical critical field Hd c - Evidently, it is connected with a symmetry 
violation of a state by the bias current (3. In Fig. |3] two stable one-fluxon states at (3 — and 
/? = 0.1 (L = 5, Hq — 1.4) are shown. It is seen that the state with (3 = 0.1 is asymmetric. 
The dynamical critical field at L = 5 are Hdc = 0.67 at (3 = and Hdc = 1-4 at (3 = 0.1. 
Upon increasing L the value of Hd c is changed {(3 = 0): Hd c (L = 5) = 0.66, Hd c (6) = 0.4, 
HdciJ) = 0.26, Hdc(8) = 0.15, if<fc(10) = 0.06, i.e. the if^c decreases. In this case the critical 
field # c i has the values: H cl (L = 5) = 1.57, tf c i(6) = 1.45, iZ" cl (7) = 1.38, H cl (S) = 1.34, 
H cl (10) = 1.28, i.e., the H c i decreases also approaching to the value of H c i(oo) ~ 1.274. 



III. TRANSITIONS BETWEEN STATES 

As it has been shown in the previous section, every stationary state of LJJ, i.e., the 
solution of Eqs. (|l])-(0), corresponds to a minimum of the thermodynamic Gibbs potential 
and these minima are not equivalent with respect to the problem of instability. For example, 
in Fig. | stationary states of LJJ at H = 2.035, (3 = 0.001 and L=5 are shown. The values of 
the Gibbs potential calculated using Eq. (f|) are as follows: Gi = —5.03, G 2 = —4.52, G 3 = 
— 4.61,(j4 = —4.64, G 5 = — 4.61,G 6 = —6.7. States 4 (Meissner) and 6 (one-fluxon) are 
stable, the other ones are metastable. It should be noted that unstable state 1 corresponds 
to deeper minimum than the stable state 4. This property contradicts the naive idea that 
more stable states occur at deeper minima. Now we shall consider this question in detail. 

The sine-Gordon equation with dissipation and bias current describing an evolution of 
initial state has the form: 

<Ptt{x, t) + ^(p t (x, t) - <f xx (x, t) = - sin <f(x, t) + (3, (7) 

where t is a time normalized to the inverse of the Josephson plasma frequency ujj = 
y / 27rcj c /C < l ) 0) C is the junction capacitance per unit area, 7 = ^qujj / AncRjc is the dis- 
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sipative coefficient per unit area, R is the resistance of junction per unit area. We write 
down the boundary conditions for Eq. ([!]) in the form 

<p x {x, t) \ x=0 = H(0, t) = <p x (x, t) \ X=L = H(L, t) = H (l- ae~ t/2t ° cos 0.5t) . (8) 

The integration of Eqs. (0)-@ for H = 2.035, (3 = 0.001, L = 5 (the same as in Fig. 
D and 7 = 0.26 gives: the metastable state 1 passes to the stable state 6 at any values 
of perturbation parameter a, 2 — > 4 at a=0, 2 — > 6 at a=l, 3 — > 4 at a=0.05, 3 — > 6 at 
a=0.07, 4 — > 6 at a=0.5 and so on. Every transition from the metastable state to the stable 
one, m — > n, is a transition from the state with the certain value of local minimum G m to 
another state with smaller value of minimum G n . These transitions m — > n with G m > G n 
are realized by certain values of the parameter of the initial perturbation a in expression 
(§). One can say that the local minima of Gi are connected with each other by a certain 
disintegration channel along the coordinate a. From this point of view one can say also that 
stationary states contain a specific "latent" parameter, by which a connection with different 
local minima Gi may be realized. In particular, the perturbation parameter a appears here 
as a "latent" parameter. It is possible, there are several "latent" parameters connecting the 
stationary states. One of the most important characteristics of "latent" parameters is that 
the stationary state does not depend on them directly; however, the form of the asymptotic 
state and the rate of disintegration depend essentially on them. The presence of a "latent" 
parameter apparently explains, a nonequivalence of the different local minima with respect 
to the stability, especially in the case when a stable local minimum is above of a nonstable 
local one. In Table | results of the integration of Eqs. flU), (0) and the calculation of G 
for every of these solutions at Ho = 1.174, (3 = 0, L = 8 are represented. The transitions 
between states m — > n are defined as follows: the mth solution of the stationary Ferrell- 
Prange equation (^)-(fJ) was taken as an initial condition of the sine-Gordon equations ([7p 
and (H). If this mth state was unstable then it fell into the nth stable state. 

The scheme of the transitions between states m — ► n is represented in Fig. [5| It is seen, 
that G m > G n for all the transitions (we note that G3 and Gg for the metastable states 3 
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and 9 are less than Gu; the state 12 is stable). The stable states - Meissner, one-fluxon, 
and two-fluxon — are shown in Fig. [] at the same parameters as in Fig. ^. 

IV. CHAOS STRIP 

As we noted above, a number of stationary states decreases with approaching to the 
bifurcation line 0-2, but the number of nonstationary asymptotic states is increased si- 
multaneously. Changing the perturbation parameter a we can obtain three sorts of typical 
states: stationary, regularly and chaotic.0 These states are distinguished not only by a form 
of the field distribution in the junction and a variation in time, but also by values of the 
Lyapunov exponent A: for the stationary states A < 0, for the regular states A < and for 
the chaos states A > 0. The Lyapunov exponents were calculated in the same way as in 
Ref. 13. However, as the calculations have shown, chaotic states may be excited not in the 
whole region 2 (see Fig. [I]), but only in the bounded region in close to the bifurcation line 
0-2. This region is extended in the form of a narrow strip along the bifurcation line 0-2 
approximately from 0.7 to 1.6 in Hq and in the range of 0.002 - 0.015 in j3. We note, that 
the chaos strip is arranged mostly under the bifurcation line in the region 2, but not in the 
region 0, as it may be expected because of all states in the region are nonstationary. The 
chaos strip is outlined on the parametric (3 — H plane in Fig. |l[ 

This chaos strip along the bifurcation line 0-2 calls to mind (to a certain extent) the 
separatrix of a nonlinear oscillator, where a chaos motion is observed. 

V. DISCUSSION AND CONCLUSIONS 

In the present work we have shown that the parametric (3 — Hq plane of a LJJ is separated 
on series of regions with the different number of solutions of the stationary Ferrell-Prange 
equation. The boundaries between these regions — bifurcation lines — characterize an essen- 
tial modification of the system. A chaos strip arises along the bifurcation line 0-2. We have 
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found that the chaos strip is arranged in the main below the bifurcation line 0-2, where 
stationary states take place. 

We have introduced the definition of a dynamical critical field as the lowest field at which 
the one-fluxon state becomes stable for the first time in the Lyapunov's sense. In addition, 
the Meissner state may also be stable at same parameters. Because both the Meissner and 
the one-fluxon states may be thermodynamically advantageous simultaneously, our definition 
based on the stability in the Lyapunov's sense characterizes an important feature of the 
stationary states of the LJJ. 

We have shown that disintegration of the metastable states and the transition to some 
stable states m — > n occur for G rn > G n . A metastable state corresponds to the local 
minimum of the Gibbs potential, and also this minimum may be lower than this one of a 
stable state. A nonequivalence of these local minima we explain by means of existing of 
a "latent" parameter not detecting in a stationary state, by which, for example, two local 
minima may be connected and a channel of the disintegration of the upper state may arised. 
In our case the perturbation parameter plays a role of the a "latent" parameter, however, 
the number of these parameters may be much greater. We note the anology between the 
quantum transitions and the transitions mentioned above, although the system is described 
by the classical Ferrell-Prange and sine-Gordon equations. 

We are aware that we could not touch upon all questions concerning the properties of a 
LJJ. We hope to return to the problems of a LJJ in our next work. 

These investigations are supported by the Russian Foundation for Basic Research (project 
N. 96-02-19321). 
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FIGURES 

FIG. 1. Bifurcation lines. The numbers of solution of the Ferrell-Prange equation (l)-(2) are 
pointed out. The number of stable states is indicated in brackets. M denotes a stable Meissner 
state and If denotes a stable one-fluxon state. L = 5. 

FIG. 2. Dependence of the magnetic field at x = L on the phase taken at the left side of 
junction x = at H = 0.5, L = 5, (3 = 0.25 and (3 = 0.45. 

FIG. 3. One-fluxon states at H = 1.4 and L = 5 for (3 = and (3 = 0.1. 

FIG. 4. Stationary states of LJJ at H = 2.035, (3 = 0.001 and L = 5. States 1, 2, 3, 5 are 
unstable, states 4 and 6 are stable. 

FIG. 5. The scheme of transitions between states m — > n. States 8, 10 and 12 are stable 
(8 — Meissner, 10 — one-fluxon, 12 — two fluxon), others are unstable. Hq = 1.174, L = 8, (3 = 0. 

FIG. 6. The stable states: M — Meissner, If — one-fluxon, and 2/ — two-fluxon at the same 
parameters as those in Fig. 5. 
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TABLES 



Number of state 



Stability 



G 



Transitions m — ► n 



Sort of stable states 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 



unstable 
unstable 
unstable 
unstable 
unstable 
unstable 
unstable 

stable 
unstable 

stable 
unstable 

stable 



2.34 
2.78 
0.64 
14.69 
14.98 
14.69 
13.53 
-1.42 
0.64 
-0.44 
2.34 
2.29 



I -► 10 
2^8 
3^8 

4 -> 10 

5 -> 10 

6 -> 10 
7^ 10 
8^8 
9^8 

10 -» 10 

II -» 10 
12 -» 12 



Meissner 



1 fluxon 



2 fluxon 



TABLE I. Transitions between states 
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